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Abstract 

The drag force of charm quarks propagating through a thermalized system of Quark 
Gluon Plasma (QGP) has been considered within the framework of both conformal/non- 



p ; 

conformal Anti de Sitter (AdS) correspondence. Newly derived Einstein Fluctuation- 
Dissipation relation has been used to calculate the heavy flavor diffusion coefficients. 
Using the drag and diffusion coefficients as inputs Langevin equation has been solved to 
study the heavy flavor suppression factor. It has been shown that within conformal AdS 
correspondence the D-meson suppression at LHC energy can be reproduced where as the 
non-conformal AdS correspondence fail to reproduce the experimental results. It suggests 
collisional loss alone within non-conformal AdS correspondence can not reproduce the 
experimental results and inclusion of radiative loss becomes important. 



1 Introduction 



The nuclear collisions at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC) energies are aimed at creating a new state of matter where the bulk prop- 
erties of the matter are governed by the light quarks and gluons. Such a state of matter 
is called Quark Gluon Plasma (QGP). The study of QGP is a field of great contemporary 
interest and the heavy flavors, mainly, charm and bottom quarks, play a vital role in such 
studies. This is because heavy quark do not constitute the bulk part of the system and 
their thermalization time scale is larger than the light quarks and gluons and hence heavy 
quarks can retain the interaction history very effectively. Therefore, the propagation of 
heavy quarks through QGP can be treated as the non-equilibrium heavy quark executing 
Brownian motion in the thermal bath of QGP and the Langevin equation can be used to 
study such a system. 



In the recent past several attempt has been made to study both heavy flavor suppres- 
sion [1, 2] and elliptic flow [3] within the framework of perturbative QCD [4, 5, 6, 7, 8, 9, 
10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. However it is pointed out that the perturbative ex- 
pansion of the charm-quark diffusion coefficient is not well convergent [20] at the temper- 
ature range attainable at RHIC and LHC collisions. Hence, non-perturbative [21, 22, 23] 
contributions are important to improve heavy quark diffusion. One possible alternative 
way to estimate the drag force is the gauge/string duality [25], namely the conjectured 
equivalence between conformal N=4 SYM gauge theory and gravitational theory in Anti 
de Sitter space-time i.e. AdS/CFT. Some attempts have been made in this direction 
to study heavy flavor suppression. Within this AdS/CFT model RHIC results has been 
reproduced well [26] whereas a few other attempts [27, 28] suggest that AdS/CFT under 
predicts recent ALICE results [29]. The non-zero value of bulk viscosity obtained form 
lattice QCD calculations [30] indicate that at the temperature range relevant of RHIC 
and LHC collisions the fluid behavior is non-conformal. Therefore, it would be interesting 
to construct a gravitational dual which captures most properties of QC. This can be done 
by breaking the conformal symmetry in the AdS space and construct AdS/QCD models 
[31, 32]. In this paper we have made an attempt to test these AdS/QCD models by 
studying D-meson suppression at LHC collision energies. 

2 Langevin equation and holography 

Consider a heavy quark of mass M and energy E passing through QGP at a temperature, 
T(« M). The heavy quark suffers random kicks leading to momentum transfer q ~ T in 
a single elastic collision with the thermal bath. Hence, it requires many collisions to change 
the heavy quark momentum significantly. The dynamics of heavy quarks propagating 
through the QGP can thus be approximated as a succession of uncorrelated momentum 
kicks which leads to a Fokker-Planck equation that can be realized from the Langevin 
equation [33, 9, 10]. 



where 70 is the drag coefficient, £ is the random force and D is the diffusion coefficient. 
To study the momentum evolution, we use the Ito discretization and Gaussian white noise 
has been implemented to generate the random momentum kicks. 

2.1 Drag and diffusion coefficients in conformal holography 

AdS/CFT in its original form, relates N = 4 SYM gauge theory on four-dimensional space 
time to the IIB string theory on AdS 5 x S* 5 background, where the conformal symmetry 
of SYM gauge theory is realized in the conformal isometry of dual metric [25]. This 
correspondence can also be generalized to finite temperature, where the space-time dual 
to N = 4 SYM plasma with temperature T is a black-hole AdS with Hawking temperature 
T [34]. The metric of AdS black hole is 




<6(W) >=D l3 5(t-t>) 
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where r h = 1/(ttT) is the horizon of the black hole. According to the standard AdS/CFT 
prescription, the energy-momentum tensor of dual theory is encoded in behaviour of metric 
near the boundary. Using this dictionary we find that metric (2) is dual to a plasma with 
conformal equation of state, e = 3p. 

By studying the dynamics of trailing string in this background [35, 24], it has been 
shown that the drag force exerted on a moving quark in a static N=4 SYM plasma is 
given by 

_ dp _ irV\ 2 p_ = 

r conf — ^ — 2 SYM ]\/[ ~ conformal P \o) 

Also by investigating the fluctuations around the classical string configuration [36, 37], 
one finds the transverse diffusion coefficient as follows 



D =nV\i*Tly M (4) 



where 7 is the Lorentz factory = 1/y/l — v 2 . 

Using (3) and (4) we find the following "modified Einstein relation" between drag and 
diffusion coefficient 



D = 2MT SY My/l ^ conformal (5) 

In terms of world-sheet temperature, T s = T / ^y, the above equation takes the following 
form 

D =2ET s T conf (6) 

this is the usual Einstein relation for a relativistic particle moving in a thermal bath with 
temperature T s . So the world-sheet temperature, T s = is the effective temperature 

for a quark moving in a static plasma [38]. 

In order to apply N = A SYM results to QCD, we use alternative scheme introduced 
in[39]. According to this proposal, one equates the energy density of QCD and SYM, 
which leads to T S ym = Tqcd/3 3 . Also by comparing string prediction for quark-antiquark 
potential with lattice gauge theory we find that, 3.5 < A < 8. Therefore in terms of QCD 
temperature and coupling we have 

rpl 

r cD n/ = a— — (7) 

Dconf = 2 4l^ QCD (8) 

where a = ^ = 2.1 ±0.5. 
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2.2 Drag and diffusion coefficients in non-conformal holography 

In the previous section we have described different ways of evaluating the drag and dif- 
fusion coefficients of SYM plasma. However SYM and QCD have different properties: 
equation of state, phase transition, symmetries, etc. In particular SYM plasma is a con- 
formal fluid with vanishing bulk viscosity. On the other hand, QGP looks like a conformal 
fluid at very high enough temperature, T » T c . So it would be interesting to construct 
a gravitational dual which captures most properties of QCD [40, 41, 42]. To do so, we 
have to break the conformal symmetry of AdS space. In [40, 41] a 5-dimensional non- 
conformal gravitational model dual to QCD is proposed, where the non-trivial profile 
of dilaton breaks down the conformal symmetry. The action of 5-dimensional Einstein- 
dilaton is given by 

S = -y^r / d"V=0 (R - ^(«90) 2 + V{4>)) (9) 

where G5 is the five-dimensional Newton constant. By choosing a suitable scalar potential 
one can mimics the QCD equation of state and other thermal properties. We choose the 
suggested potential in [43]: 

V{\) = ^|{1 + V X + V 1 \l[\n(l + V 2 \l + V 3 X 2 )]^} (10) 

with 

8 /23 + 364\ 2 22 51 

Vo = _ A , v 2 = % [-^j , A = _, A = (11 ) 

It has been shown [43] that this potential reproduces the lattice EOS and velocity of 
sound. Drag force in this model is calculated in [43]: 

U g2A s (r s ) 

Fnon—conf ~Z 7^ = ^ non—conf V iXty 

where v is the speed of quark, r s is the world-sheet horizon and A is a component of dual 
metric. Here l s is a free parameter which can be fixed by matching the string tension to 
the string tension derived from the lattice QCD calculations and is given by [43] 

l s ~ 0.15 (13) 

Note that unlike conformal case, non-conformal drag coefficient is velocity dependent 
through r s . It is useful to study the ration of drag force in non-conformal holography to 
the conformal case: 

/•' , I' , 9 1 

R(V,7F) ( 14 ) 



Fnon—conf ^non—conf 2.1 „ / T 



Fconf Tcon/ ^ T c 

where R is a function of temperature and velocity of quark. We reproduced this function 
in Fig. 2 and Fig.3 for completeness [43]. If we take l s as a free parameter, then this 
relation takes the following form : 
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T non - conf =^r 2 ^R(v,^)T conf = (°-^r 2.1 R{v^f^ (15) 
The modified Einstein relation (6) for non-conformal becomes 



D n on—conf 2£/ Ts,non—conf ^non—conf (16) 

where T Sjnon _ cori / is the world-sheet temperature in non-conformal case. In terms 
of the ratio of world-sheet temperature in non-conformal to conformal case, G(v , = 
s,non~conf ^ ^ e above equation takes the following form [38] 

s,conf 

T 2M _ T 

D n on-conf = 2-E T S:Con f G(v, — ) T non _ con f = — — ^7 Tqcd G(v, — ) T non ^ con f (17) 

Jc 3* J c 

There are several limits on AdS/CFT results discussed here. In [44] the effects of 
hydrodynamic expansion of QGP on drag force exerted on a moving quark have been 
studied. It was shown that there is an upper bound for velocity of quark (vi, oun d ~ 0.98) 
such that below this bound, drag force acting on the quark is just the localized version of 
static plasma (replacing the temperature in drag formula of static plasma by instantaneous 
temperature of QGP). On the other hand, for a fast quark with a velocity bigger than the 
above bound (for charm quark this bound in velocity corresponds to around lOGeV in 
energy) drag force is not a local function of the medium variables. Thus, for p? ^> lOGeV 
the local approximation is not valid due to hydrodynamic expansion. 

Also it has been shown in [38] that for p? > lOGeV the white noise approximation 
in Langevin equation breaks down. Above this bound, one needs the full frequency- 
dependent correlators to study the diffusion process [45]. 

3 Initial condition and space time evolution 

After obtaining the drag and diffusion coefficients from the conformal and non-conformal 
holography, we need the initial heavy quark momentum distributions for solving the 
Langevin equation. In the present work, the Pt distribution of charm quarks in pp colli- 
sions have been generated using the POWHEG [46] code, implementing pQCD at NLO. 
With this initial heavy quarks momentum distribution, the Langevin equation has been 
solved . We convolute the solution with the fragmentation functions of the heavy quarks 
to obtain the pr distribution of D mesons. For heavy quark fragmentation, we use the 
Peterson function [47]. 

The experimental interest is the nuclear suppression factor (Raa), defined as 

dN Au+Au 

= lff^j*i < 18 > 

JVcoll X d 2 pTdy 

a ratio that summarizes the deviation from what would be obtained if the nucleus-nucleus 
collision is an incoherent superposition of nucleon-nucleon collisions. In Eq. 18 iV co ii stands 
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for the number of nucleon-nucleon interactions in a nucleus-nucleus collision. In the 
present scenario the variation of temperature with time is governed by the the equation 
of state (EOS) or velocity of sound of the thermalized system undergoing hydrodynamic 
expansion. H 




Figure 1: Variation of (Raa) as a function of momentum px for D meson at ALICE 
within conformal model. Experimental data taken from [29]. Although we have present 
the results here upto px ~ 15 GeV, it should be mentioned here that the white noise 
approximation in the Langevin equation is not valid beyond pt = 10 GeV. In this backdrop 
the theoretical results should be taken. 



The system formed in nuclear collisions at relativistic energies evolves dynamically 
from the initial QGP state at temperature, Tj to the quark-hadron transition temperature, 
T c . Boost invariance Bjorken [48] model has been used for the space time description of 
the QGP. It is expected that the system formed in nuclear collisions at LHC energy in 
the central rapidity region is almost free from net baryon density. Therefore, the equation 
governing the conservation of net baryon number need not be considered in the present 
case. The total amount of energy dissipated in the system by the charm quarks depends 
on the number of interaction it undergoes i.e. on the path length (L) it traverses within 
the medium. The value of L in turn depends on the spatial co-ordinates, (r, 0) of the point 
of creation of the charm quark. The probability, P(r, 0) that a charm quark is created at 
(r, 0) depends on the number of binary collisions at that point. P(r, 0) is given by: 

P{r^) = ^{l-^)B{R-r) (19) 

where R is the nuclear radius. In obtaining the above expression for P(r, 0) spherical 
geometry has been assumed, therefore, it is better applicable for central collisions. The 
charm quark created at (r, 0) in the transverse plane of the medium will propagate a 
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Figure 2: Variation of R(y,T/T c ) as a function of momentum p T . 



path length, L given by L = y/R 2 — r' 2 siri 2 <f> — rcoscf). The geometric averaging has been 
performed for the the drag and diffusion coefficients along the path length. The initial 
temperature (Tj) and thermalization time (r^) of the background QGP are constrained by 
the following equation: 

3 2tt 4 1 1 dN 

1 Ti ~ 45C(3) 4a eff ttR 2 ~~dy ' ( ' 

where (dN/dy) is the measured hadronic multiplicity, £(3) is the Riemann zeta function 
and a e ff = n 2 g ef f/90 where g eff (= 2 x 8 + 7 x 2 x 2 x 3 x N F /8) is the degeneracy of 
quarks and gluons in QGP, A^=number of flavors. Initial condition for LHC energy has 
been taken from Ref [49]. 



4 Results 

With the formalism discussed above the result for (Raa) are shown in Fig 1 for the 
conformal holography. It is found that the ALICE data can be explained reasonably well 
for a = 2. For a = 3 we under predicts the experimental data. Here it may be mentioned 
that within the conformal holographic model the RHIC results was explained reasonable 
well for a = 2 — 3 [26] in Langevin dynamics. However the conformal holography model 
under predict the recent ALICE data [29] presented in ref [27, 28]. As the conformal results 
are aways from reality, in a very first attempt we are implemented the con-conformal 
results with the Langvine equation to study the D-meson suppression at LHC energy. 
The ratio of interaction obtained from non-conformal to conformal is displayed in Fig 2 
with respect to p F - It is observed that the non-conformal drag force reduced by a factor 
3-4 time than the conformal case. The momentum dependent is also weak as showm in 
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Figure 3: Variation of R(v,T/T c ) as a function of temperature T. 



Fig 2. Similarly in Fig 3 the variation of the ratio from non-conformal to conformal is 
plotted with respect to T for fixed pt- It is found that non-conformal drag force reduced 
by a factor 3-4 time than the conformal. In Fig 4 the variation of (Raa) has been shown as 
a function of Pt for various value of l s within the non-conformal holography . It is found 
that non-conformal drag force over predicts the data for a realistic value of l s = 0.15. It 
is quite expected as the non-conformal drag forces suppressed by a factor 3-4 as compare 
to conformal case (in Fig 2 and Fig 3). Since we are only considering the collisional loss 
within the non-conformal holographic it fail to reproduced the experimental data. The 
results will improve if the radiative loss from the non-conformal holography will be taken 
into account. 



5 Summary and conclusions 

In an attempt, we have studied the D- meson suppression at LAC energy within both 
the conformal and the non-conformal holographic model. It is observed that the non- 
conformal holographic model over predicts the ALICE data where as the data can ex- 
plained reasonably well for a = 2 within the conformal holography. This is because the 
non-conformal drag force suppressed by a factor of 3-4 compared to the conformal case. 
It is expected that inclusion of the radiative loss from the non-conformal side will improve 
the results. Therefore, more systematic studies are needed from the non-conformal side 
like inclusion of radiative loss, etc to improve the description of the experimental results. 
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Figure 4: Variation of (Raa) as a function of momentum for D meson at ALICE within 
non-conformal model. Experimental data taken from [29]. 
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